Why Noether symmetry of F(R) theory yields three-half power law? 
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Abstract 



Noether symmetry of F(R) theory of gravity reveals F(R) — , if the expression for the scalar curvature 
for R-W metric is treated as a constraint and entered into the action through a Lagrange multiplier. In the 
process, a cyclic coordinate is found which gives a solution that appears to explain the present cosmological 
evolution. The interesting issue is that out of infinite curvature invariant terms, Noether symmetry selects only 
R? . Here, we explore the very speciality and study the attractive features and shortcomings of such term in 



C/3 , 

the context of cosmological evolution. 
PACS 04.50. +h 

> : 

1 Introduction 

00 . 

, In the last decade several path finder models have been explored which enunciated the best fit with SNIa data set. 
For most of these models SNIa data set implies a smooth transition from a decelerated phase to a recent accelerated 
phase of cosmic evolution, and the presence of nearly 74% of dark energy, in disguise. However, starting from the 
^vq ' ACDM concordance model, often cited as the standard model, all the models suffer from some sort of disease. 
t-H . The order of magnitudes of the vacuum energy density (corresponding to cosmological constant A ) calculated by 
the Quantum field theorists and the Cosmologists differ by a huge (10 120 ) amount. Quintessence model solves 
the problem, but it requires a scalar field with a suitable potential. But then, we need a scalar with some exotic 
potential, which is not at hand presently. Further, Quintessence model does not suffice to handle the situation 
with state parameter uj < — 1 , which is not excluded in view of the presently available data. It requires additional 
scalar field containing kinetic term with reverse sign. Thus, the situation gets more and more complicated. There 
was attempt to solve the puzzle assuming particle creation phenomena at the expense of gravitational field [TJ , [5] , 
which has strong theoretical basis, but there is no direct experimental evidence as yet. The only thing that is left 
is curvature, which can not be detected directly. Inflation in the early Universe may be an artefact of R 2 term in 
gravity [3] and comic acceleration in the late Universe might have been caused by term [4]. Thus, it appears 
that an action in the form S = J [ctR + PR 2 + ^R^ 1 ]^ — gd 4 x , having pure gravitational origin, might be the one 
to solve the cosmic puzzle from early Universe till date. 

However, term fails to produce Newtonian gravity in the weak energy limit and does not comply with the 
solar test [5]. Even worse is, it shows unavoidable instabilities within matter in the weak gravity limit [6] and 
fails to explain big bang nucleosynthesis (BBN) [7J. Further, for any R~ n term, with n > 0, it has been shown 
[8] that the scale factor during the matter phase grows as 1 2 , instead of the standard 1 3 - law, prior to the 
transition to an accelerated expansion. This is grossly inconsistent with cosmological observations viz., WMAP 
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(Wilkinson Microwave Anisotropy Probe) and LSS (Large Scale Structure). Therefore one requires to find some 
scalar invariant other than , suitable to explain late time cosmic acceleration and that may pass the solar 
test. Such a theory, dubbed as F(R) theory of gravity [see [9] and references therein] is much more attractive, 
promising and intriguing than scalar held models. But then, there are at least a couple of associated problems. 
Firstly, out of indefinitely many curvature invariant terms, one requires an elegant method to choose suitable ones 
that fit all cosmological data at hand. Next, it is practically impossible to solve the field equations of F(R) theory 
of gravity, corresponding to any nonlinear form of R . Both the problems presumably may be solved by invoking 
Noether symmetry. 

In order to express higher order theory of gravity in non-degenerate and canonical form, it is required to be 
spanned by an auxiliary variable in addition to the basic variables. In the process, the configuration space of the 
Lagrangian gets enlarged. Thus, one can demand the existence of Noether symmetry amongst the field variables 
(viz., between the scale factor and the auxiliary variable) even when the theory is reduced to a measure zero 
subset of its configuration space, i.e., if the Robertson- Walker minisuperspace model is accounted for. This was 
attempted largely by several authors [10], [11], to find a suitable form of F(R), spanning the Lagrangian by a 
set of configuration space variables (a, a, R, R) . Noether symmetry of F(R) theory of gravity [10], [11], yields 
F(R) = i?2 , in the Robertson- Walker minisuperspace model, which appears to be important in studying early 
Universe but not the late. However, explicit solution indicates a smooth transition from early deceleration to the 
late time cosmic acceleration [11] , which is definitely an attractive feature of such a curvature invariant term. 
Nevertheless, the more interesting issue is that, out of infinite number of curvature invariant terms, Noether 
symmetry selects R? term in particular. The result is doubtful for the following reason. In a series of works 
|12j . [13] it has been pointed out that although classical field equations can be formulated choosing any auxiliary 
variable, different auxiliary variables produce different Hamiltonian. As a result, different quantum descriptions are 
revealed. A toy model [13] explains the reason, which is, different auxiliary variables produce different definitions of 
canonical momenta. In the classical field equations, it does not make any difference, since one requires derivative of 
momenta. But Hamiltonian depends on momenta itself, and is different for different choice of canonical momenta. 
In finding Noether symmetry, the situation is the same. Here again we require momenta, and not its derivative 
and so different auxiliary variable might produce different results for Noether symmetry. The standard technique 
of finding auxiliary variable is to take partial derivative of the action with respect to the variable with highest 
derivative appearing in the action, i.e., Q = -QM- for the Robertson- Walker line clement, where hij is the metric 

on the three space. Thus, we think that the study of Noether symmetry of F(R) theory of gravity requires further 
attention. 

Despite the fact that i?2 theory of gravity has some attractive features, particularly in explaining the late time 
cosmic evolution it suffers from the same type of problems associated with R~ n , with n > theories of gravity, 
since it never ends up with Newtonian analogue in the weak field limit. This is also apparent from the solution 
obtained |llj . The solution depicts that in the matter dominated era, prior to the transition to an accelerated 
expansion, the scale factor never follows the standard t~s -law. Such disease may be cured if such an action would 
have been supplemented by Einstein-Hilbert action. Hence, in the following section we try to explore Noether 
symmetry of F(R) gravity being supplemented by the Einstein-Hilbert term. This is done to understand if 
the technique of finding such symmetry is correct. If it is, it is then supposed to yield the same result, since 
aR + j3F(R) = jFi(R) , and this is what we have obtained. For a further check, in section 3, we explore Noether 
symmetry of i?2 term following the standard technique and end up with the same solutions obtained earlier [TT]. 
Thus we understand that Noether symmetry of pure gravity, yields nothing more than F(R) = R^ . This initiates 
to find what is so special in such a curvature invariant term. In section 4, under a change of variable hij , we 
found that Noether symmetry is in-built in Ri term, since hij = a 2 becomes cyclic and the field equations are 
solved at ease yielding the same known solution. This implies that one does not require to request [11] Noether 
symmetry of R? curvature invariant term. In section 5, we review Noether symmetry of F(R) theory of gravity, 
in view of the changed variable hij , and observe that the conserved current may be solved quite easily and it is 
not required to find cyclic variable any more. In the next section 6, we introduce matter and found solutions both 
in the radiation and in the matter dominated era. Section 7 is devoted to explore how far R^ term is acceptable 
in the context of available cosmological data. In section 8, we conclude with our findings. 



n 



2 In search of Noether symmetry for Einstein-Hilbert action being 
modified by F(R) theory of gravity. 



We have already discussed the problem of F(R) theory of gravity if it does not contain Einstein-Hilbert part in 
the action. Noether symmetry of F(R) theory of gravity yields F(R) = R? . The solution is encouraging as far 

2 

as late time cosmic acceleration is concerned. But the problem is that, it only passes transiently through a cx 1 3 
in the matter era before the transition to accelerated phase [TT], and thus fails to explain WMAP and LSS, as in 
the case of R~ n term, [7], [8]. So, we try to find Noether symmetry for the F(R) theory being coupled to the 
Einstein-Hilbert sector. Although it appears to be trivial, since these two terms together may be combined to get 
yet another F(R) theory, however, it is important to see if it is true while Noether symmetry is being explored. 
In the process, it turns out to be an important test to check if the technique of exploring Noether symmetry is 
correct. Hence we start with the following action, 
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(1) 



where, ^ 



d 3 x, is the surface term 



§[^fhK]d 3 x is the Gibbons-Hawking boundary term and £ = 2B J VhF jR K 

for F(R) theory of gravity that emerges under variational principle (metric formalism). In the above, h is the 
determinant of the metric on three space, K is the trace of the extrinsic curvature and F tR is the derivative of 
F(R) with respect to R. Note that this surface term £ reduces to the Gibbons- Hawking boundary term taking 
F (R) = 16 \ G R and that for the curvature squared action F(R) = R 2 [15]. Action (1) in the background of 
isotropic and homogeneous Robertson- Walker line element, 
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a 3 dt d 3 x 



8nG 



[VhK]d 3 



2B 



VhF R K 



d 3 3 



(3) 



where, the expression for R has been introduced as a constraint through Lagrange multiplier A. Now varying the 
above action with respect to R, one gets A = BF, R , which is substituted in the action. Now under integration 
by parts the surface terms get cancelled and the action may be expressed as, (taking = M), 



A 



-Mad 2 + Mka + Ba 3 (F - RF, R ) - 6Bad 2 F, R -6Ba 2 aF, RR R + 6BkaF, R 



dt. 



and the Noether equation may be written as, 

XL = £ X L = a(-Md 2 + Mk + 3Ba 2 (F - RF, R ) - 6Bd 2 F, R -12BadRF, RR +6BkF, R 
+(3{~Ba 3 RF, RR -6Bad 2 F lRR -6Ba 2 aRF, RRR +6BkaF, RR ) 

'd/3. d/3 



-^d + -^R ) (-2Mad - 12BaaF, R -6Ba 2 RF, RR 
oa oR 



(4) 



»)a a+ 8R R ' (- 6Ba2hF >R* ) = °> ( 5 ) 



where, X is the vector field. Now we equate the coefficients of a 2 , R 2 , dR and others respectively to zero and 
obtain, 



(M + 6BF, R )(a + 2aa ia ) + 6B(/3 + a/3, a )aF, RR = 0. 



(6) 



QBa 2 F , RR a, R = 0. 



(7) 



6B/3a 2 F, RRR +6B(2a + aa, a +a(3, R )aF, RR +2(M + 6BF, R )aa, R = 0. (8) 



a[Mk + 3B(F - RF, R )a 2 + 6BkF, R ] + B[3[Qk - a 2 R]aF, RR = 0. (9) 

In the above equations comma ( , ) represents derivative. Equation (7) implies a ^ a(R) , since, we require 
nonlinear form of F(R), ie., F lRR ^ 0. So, we have, a = a(a) , and separating (3 = A(a)S(R). Equation (8) is 
expressed as, 



F, RR + dR~ A 



a da 
a da 



= ci, (10) 
where, C\ is a separation constant. The above equation (10) yields, 



da „ a , 

+ 2- = - Cl A, 11 
da a 

and, 

F,rrr S, R Ci 



F, RR S S ' 
Equation (6) can as well be separated using the separation constant C2 as, 



(12) 
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S= 6B Z in , (13) 



iR 
C2F, RR 

and 



a 2 A, a +aA + c 2 (2aa, a +a) = 0. (14) 
Equations (11) and (14) are combined to give the following differential equation, 

a 2 a, aa +(3 - 2c 1 c 2 )aa, a -c x c 2 a = (15) 
In view of equations (13) and (15), equation (9) takes the following form, 

( AI + QBF ) 

k(M + 6BF, R )[( Cl c 2 - 2)a - aa, a )] + 3B Cl c 2 (F - RF, R )aa 2 + ^ ^(aa, a +2a)a 2 R = 0. (16) 

Thus, one can solve a in view of equation (15) and hence A from (11). F(R) and S(R) may be solved in view 
of equations (12) and (13). Finally, equation (16) is satisfied for cic 2 = 2. At the end, following set of solutions 
exists for k = ±1, 0. 

n „ 2nR _ 1 , • a 3 d ( R 



where, n is a constant of integration and F is the conserved current. However, the most interesting result is, 

F{R)^-^R + CRi 1 (18) 



which, under the substitution M 



SttG 



, reduces the action to 



-g d 4 x[DR^} 



(19) 



Thus, Noether symmetry of F(R) does not admit an explicit Einstein-Hilbert term in the F(R) theory of gravity. 
This is the result that we may expect if the technique of exploring Noether symmetry, treating the expression of 
R as a constraint is correct. So, in the following we further check, what Noether symmetry yields, starting from 
i?2 term, instead of F(R). 
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3 Noether symmetry of - standard technique. 

This section is important for the following reason. To express F(R) theory of gravity in non-degenerate and 
canonical form, one usually introduces the form of R as a constraint through a Lagrange multiplier A. The 
extended configuration space in such case is (a, d, R, R) . In the standard formalism, where a form of curvature 
invariant term is fixed, the configuration space is different viz., (a, d, Q, Q) . Here, Q is an auxiliary variable 
chosen as the derivative of the action with respect to the highest derivative of the variable appearing in the action, 
that has not been possible to integrate by parts to produce a surface term. We here check, if Noether symmetry 
along with the solutions presented by [H] exists, when action containing i?2 term is spanned by such a set of 
configuration space variables, i.e., whether Noether symmetry is an artifact of a particular choice of configuration 
space variables or not. One may expect that since Noether symmetry of F(R) theory yields i?2 , so £xL = 
will be exactly solved if F(R) = R^ is assumed a-priori. However, we show below that this is not the situation. 
The action under consideration in the Robertson- Walker minisuperspace model is, 



S = B Ri^/^jd 4 x -2B \VhF, R K~\ d 3 x = B / 6^ [ad + d 2 + k]%dt - 9V6B[aa + d 2 + k\^aa 



Now let us define the auxiliary variable Q as, 
dS 

Q = — = 9V6Bfad + d 2 + fcl^a. 
oa 

Expressing the action in non-degenerate, canonical form 

l Q 3 
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and integrating by parts, we obtain, 
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(21) 



(22) 



(23) 



One can easily check that the last two terms get cancelled and so, there is no controversy on the choice of the 
auxiliary variable. The field equations are, 
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1458S 2 a 3 

Equation (24) is a fourth order equation in the scale factor (a), while equation (25) is third order and apparently 
it is impossible to solve the pair. So, it is better to ask Noether symmetry between the field variables a and Q, 
i.e., £xL = XL = 0, which reads, 



486B 2 a 4 



-Q 
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P a a + /3 Q Q)a = 0. (26) 



Equating coefficients of a 2 , Q 2 , aQ and terms free from derivative, we get following four equations, which are 
linear in a and j3 , viz., 



aQ — 2aa t aQ — a/3 



a 2 /3, a 



(27) 



= 0. 
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In the above comma represents derivative, as before. Equation (28) dictates that a ^ a(Q) , while equation (30) 
clearly yields two cases, viz., 



-Q, and, Q = 9\/6BVfc 



3.1 Case - I, p = f Q 

Equation (29) in view of the above condition yields,, 



a = — , and so, (3 — —^Q- 



(31) 



It is quite easy to verify that equation (27) is trivially satisfied for the above solutions. Thus the conserved current 
T is, 



T = ap a + j3pn =>■ — [ — ] = Constant. (32) 
dt \a J 

Plugging in the definition of Q from relation (21), one can immediately solve the above equation for a as, 



a 4 i + a 3 t + 



3a| 

8d4 



k I t + a2t + ai 



(33) 



where, a\ , a 2 , a% and a^ are integration constants. It is not difficult to check that the above solution satisfies 
both the field equations (24) and (25). This is exactly the solution obtained and analyzed by Cappozziello-et-al 
[11] , Note that, the solution has been found directly from Noether symmetry as we found earlier in the case of R 2 
theory of gravity being coupled with a scalar field [14] and Gauss-Bonnet dilatonic interaction [16] and we don't 
require to solve the field equations any more. This is a unique feature of higher order theory of gravity. The above 
vacuum solution may be treated as one corresponding to the early Universe, when curvature played the dominant 
role. The Universe started evolving as a oc \ft, and soon transits to a power law inflationary era. 

3.2 Case - II, Q = 9V6kBa 

Since Q is proportional to a, so, the field variables practically reduce to one, and it is nonsense to try to find 
Noether symmetry. Rather, we can check, if such a relation between Q and a yields a solution in view of the 
definition of Q given in (21) that satisfies the field equations (24) and (25). The solution is, 



\J dit + d , 



(34) 



where, di and do are integration constants. This solution satisfies the field equation (24), but equation (25) is 
satisfied only under the condition, k = . This is trivial, since then Q — and the field equations are trivially 



satisfied. In view of the definition of Q, the Ricci scalar R — 0, and so, it automatically leads to Friedmann 
solution in the radiation dominated era in the following manner. The combination of Friedmann equations, 



(a 2 k 

2- + + = -%irGp and 3 — + — I = SnGp, 
a a z a z \ a z a z 



together yield, 

+ ^ + =R = 8irG( P - 3p). 
\a a z a z J 

Now, in the radiation dominated era, p — 3p = 0, and so the above equation reduces to, 

„ fa a 2 k \ „ 
6(- + -=. + -5. =i2 = 0. 

\a a A a z / 

Further if k = 0, the solution of the above equation is a oc fa . 



4 What is so special in R^l A different look with change of variable. 

So, we observe that the technique of exploring Noether symmetry for F{R) theory of gravity treating the expression 
of R as a constraint of the theory appears to be alright, since standard technique does not yield anything else. 
Hence, the question that naturally arises, what is so special in the curvature invariant term i?s . In this section, 
we pose to answer this question. It has already been emphasized that the basic variable is hij and not the scale 
factor a . In the case of higher order theory of gravity, the configuration space gets enlarged and the new variables 
that finally come into rescue is [13] . So, here we start with the basic variable, hij — a 2 — z. In view of the 
expression of Ricci scalar, 



R = 6 . 

2z z 
The action now reads 



S = B J 3\/3(z + 2k)^dt - 2B J VhF M K 
The auxiliary variable is 



d A x. 



(35) 



^ dS 9^3 „ n . 
az 2 

and the canonical form of the action is 
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The definition of the auxiliary variable is restored immediately from the Q variation equation, while the variable 
z turns out to be cyclic. Thus the field equations are oversimplified to, 



= o, 
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zQ + 2kQ =0, 
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In view of the definition of the auxiliary variable Q given in (36), equation (38) immediately gives the following 
solution, 



a 4 r + a 3 r + 



^ k ) / ; ' • 

8a4 



(40) 



which resembles exactly with the solution (33) obtained earlier in view of Nocthcr symmetry. Note that, equation 
(39) also yields the same above solution. Thus, it is clear that Noether symmetry is in built in i?5 theory of 
gravity, since z turns out to be cyclic automatically. 

5 Revisiting Noether symmetry of F(R) theory of gravity 

We now understand that to find Noether symmetry for F(R) theory of gravity, it is better to choose the variable, 
hij = a 2 = z , rather than the scale factor a . In this section we show that under such a choice of variable, Noether 
symmetry of F(R) theory of gravity yields a rather handy conserved current, which can be solved at ease, and it 
is not required to search for the cyclic co-ordinate. 
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(41) 



Following the same technique, i.e., varying the action with respect to R and setting it to zero, one finds A = F iR , 
which when substituted in (41) and the surface term is accounted for, through integration by parts, the action 
may be expressed in the following non-degenerate and canonical form, 
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Demanding Noether symmetry one obtains, 
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Equating coefficients as usual, following four equations emerge, 
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For a non linear form of F(R), i.e., F tRR ^ 0, equation (45) implies that a — a(z) only. Thus, equation (44) 
demands separation of variable as, /3 = (3i(z)(32{R) , and is expressed as, 



while, equation (46) is expressed as, 
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C2, 



,RR 



(49) 



where, c\ and C2 are separation constants. It is not difficult to check that middle terms of the above two equations 
arc the same, i.e., @2 F p iR — P2 F p^ R + &2,r and so, c\ = C2 = c. Thus equations (48) and (49) may be solved 
and one gets, 



a = mz + n and /3i 



1 

2cz 



(3mz + n). 



In view of the solutions (50), equation (47), after simplification takes the following form, 



3(mz + n)F - 12kmF jR - 2nRF jR , 



(50) 



(51) 



which is satisfied only under the condition m = 0. Hence one can solve equation (51) for F(R) and then equation 
(48) for 02 (R) ■ The final set of solutions thus obtained is, 
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a = n, (3 = , F(R) = R*, 

z 



(52) 



where, n is a constant. The conserved current now reads, 



_ 9n d , 1 — . 



(53) 



which is having the same form as obtained earlier by Capozziello-et al [llj . Substituting the expression for R, the 
above conserved current can now be solved directly and one obtains, 



'3a 2 

a 4 t A + a 3 t 3 + ( ^--k)t 2 + a 2 t + a 1 , (54) 
which is the same solution obtained earlier in (33) and (40). 

6 Presence of matter 

The solution (33), or ( (40) or (54)) , that we have presented so far corresponds to vacuum, which may have some 
importance in the early curvature dominated Universe, since it represents transition from decelerating Universe 
with a oc y/t to a power law inflation. To unveil the importance of R? in the late stage of cosmic evolution, one 
has to take some form of matter into account. So, let us start again from the action, 
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BRi - Mo _3( ™ +l! 



V5d 4 x - 2B 



-VhR K 
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d A x, 



(55) 



where, w — - is the state parameter. Note that, in the radiation dominated era, w = ^ , and so, the second term 

reads, Ma~ 3 (' w+1 ^ = and in the matter era, w = 0, and it reads Mo" 3 ''° +1 ' = ^r- ■ Here, p r o, and p m o 
refer to radiation density and the matter density at the present epoch, which has together been represented by 
M . Following the same prescription of introducing auxiliary variable given in (36), and after taking care of the 
surface term, the above action can be expressed in the following canonical and non-degenerate form, 



S = 



-Qz + 2kQ 



729B 2 



Q 3 - Mz~ — 
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(56) 



The field equations are, 



Q + — -z-^— >=Q, (57) 



Qz + 2kQ-^^Q 3 = Mz-^, (58) 
In the radiation dominated era w = 3 , and the field equations admit a solution in the form, 

= a = 6 ^ + &i, (59) 

under the condition, fc = and £ = — , instead of the standard Friedmann type, i.e., a oc t% . This creates 

problem in explaining BBN, since the Universe expands at a much faster rate than the standard Friedmann model. 
Further, due to such a high expansion rate in the radiation era, it becomes difficult for the seed of perturbations 
to accumulate matter, in the matter dominated era, to form structures. Now, in the matter dominated era w = , 
and so is Q . Thus we recover the same solution (33) or, ((40) or (54)) , i.e., 



3a 2 

a 4 t 4 + a 3 t 3 + [ g- 3 - - k ) t 2 + a 2 t + ai 



(60) 



which is true for arbitrary curvature parameter, viz., fc = 0, ±1. 



7 Fitting the observational data 

Data fitting of the above unique solution for Ri in the matter dominated era, has been presented already by 
Capozziello-et al [11]. Setting H to = 1 and deceleration parameter q = —0.04, they reduced the solution to 
a one parameter model, for which a nice fit with ACDM is obtained. However, for such parametric values, the 
rcdshift at decoupling, Zdec — 1100, occurs at t = 11 X 10 6 years, which is clearly in contradiction with the WMAP 
data and the standard model. Further, the Hubble parameter H(z) versus redshift z curve (figure 1), is largely 
different from the standard model at high redshift. For the required fit, they [IT] also took some finite value of fc, 
viz., k — —0.49. Here we observe that the only power law solution that satisfies the field equation in the radiation 
dominated era, requires k — 0. So, we are required to set k = 0, in the matter dominated era also. In the solution 
above, one can set a\ = 1, without any loss of generality. Still, we are left with three parameters, viz., a 2 , a 3 and 
a 4 . We have studied following two cases. 

7.1 #0*0 = 1- 

Fixing up i?o^o = 1 , one can now express, one of the parameters of the theory (we have chosen a 2 ) in terms of 
the other two (viz., 0,3 and 04 here). The important observation is that the fit with SNIa data is almost perfect 
for a truly wide range of parametric values of a% and 0,4. We briefly present our result in table 1. The table 
depicts that the present value of the state parameter w e o , remains close to — 1 for a huge range of parametric 
values of 10 -3 < a 4 < 10 6 (a 4 ^ 0) and 10~ 5 < a 3 < 10 6 . Within such a range of values the transition redshift 
z a always remains close to 1 . The Redshift-Distance modulus curve fits almost perfectly with the ACDM model 
and thus the SNIa-data in disguise, taking the density parameter of the baryonic and nonbaryonic dark matter 
n m = 0.26 and the rest, = 0.74, as the density parameter of the dark energy. One of such fits has been shown 
in figure 2. Despite such amazingly attractive result, some serious problems have also been encountered. Almost 
in all the cases the redshift value is z « 200 at around t = 3.75 x 10~ 4 Billion years. Note that according to the 
standard model it is the time after Big-Bang, when photons decoupled and whose redshift value is z& ec k, 1100, 
as confirmed by WMAP data. This is practically a huge problem. Further, considerably large deviation from the 
standard model in the early matter dominated epoch is clearly visible from figure 3, where a plot of the Hubble 
parameter H(z) versus the redshift parameter z is drawn with the parametric values 04 = 1,03 = 0.1. The 
two, deviates even further at a large redshift, and the feature is the same for all other parametric values. Such a 
deviation from the standard model in the early matter dominated era is inconsistent with the WMAP data and 
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Figure 1: Hubble parameter H(z) versus the Redshift parameter z for the present model (blue, lower curve) 
shows a huge deviation at large redshift, with the standard model a oc (red) for the parameters chosen by 
Capozziello-et al [llj . 
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Figure 2: Distance modulus (M — m) versus redshift z plot of the present model (blue), shows perfect fit with the 
ACDM model (red) taking 04 = 1 and 03 =0.1. The same feature is observed almost for all parametric values of 
03 and 04 . 



LSS [5] . Finally, Figure 4 is the effective state parameter w e versus redshift z plot for the same above parametric 
values. It shows that the effective state parameter w e remains close to 0.33 at large z, which clearly depicts 
that the early matter dominated era is far from the standard Friedmann model (w e =0) and tracks a oc 1 3 only 
transiently. 

7.2 H t = 1 and z dec = 1100, for t = 3.75 x 10 5 years. 

In the preceding section we have noticed that the unique solution obtained for i?l action in the matter dominated 
era shows excellent fit with ACDM model and hence SNIa data in disguise, for a huge range of parametric values, 
but it fails to fit WMAP data corresponding to the early matter dominated epoch. Hence, here we fix up the 
WMAP data a-priori, which fixes up yet another parameter of the theory, viz., 03, in the present case. Thus, we 
end up with a one-parameter model. Here, we observe that again for a wide range of parameter, 04 > 10 -5 , with 
04 7^ 0, the Redshift-Distance modulus curve fit is almost perfect, like the one presented in figure 2. The present 
value of the effective state parameter is w e o rj —0.35 in almost all the cases, which is no where near ACDM 
model. But then, w e o is model dependent, and so it creates no problem as such. However, the transition redshift 
is unacceptably large, the minimum of which is z a = 4.2 . Further, the Hubble parameter H{z) versus redshift z 
plot (figure 5), here again shows huge deviation which increases rapidly with z. The only difference is the Hubble 
parameter for the present model here, remains far below than that corresponding to the standard model. The 
effective state parameter again tracks 0.33, instead of w e — 0, in the early matter dominated epoch, which is 
grossly inconsistent with the LSS. 
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Figure 3: Hubble parameter H(z) versus the Redshift parameter z for the present model (blue) shows a wide 
deviation at large redshift, with the standard model a oc ts (red) for Hotg = l,h = 0.72, 04 = 1 and 03 = 0.1. 
Blue being on the top. 
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Figure 4: The effective state parameter w e versus versus the Redshift parameter z for the present model (04 = 
1, as — 0.1) shows that in the early epoch of matter domination the state parameter tracks the value w e = 0.33, 
which means the model always tracks a oc 1 2 rather than the standard model a oc 1 3 . This creates problem to 
explain large scale structure. 



Table 1: Data fitting with H Q t Q = 1 and h = 0.72 
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Figure 5: Hubble parameter H(z) versus the Redshift parameter z for the present model (blue) again shows a 
wide deviation with the standard model a oc (red), while blue now is below the red. The plot is for 04 = 1. 

8 Summary 

We have made a detailed analysis corresponding to the cosmological evolution with F(R) = R% theory of gravity, 
in the Robertson- Walker minisuperspace model. We enlist the present findings. 

1. Noether symmetry of F(R) theory of gravity does not admit anything else than F(R) = R? in the Robertson- 
Walker minisuperspace model. 

2. The speciality of such a curvature invariant term is apparent under a change of variable from the scale factor 
a to the basic variable hy = a 2 , since hy becomes cyclic. Thus we conclude that Noether symmetry is in-built 
in Ri theory of gravity, for Robertson- Walker minisuperspace model. Under the change of variable, the field 
equation looks pretty simple and is solved at once. 

3. If one tries to find Noether symmetry of F(R) theory of gravity, using the basic variable h+j , the conserved 
current is solved directly, and it is no more required to find cyclic co-ordinate. 

4. The solution in the radiation era (a oc t^), is substantially different from the standard Friedmann solution 
a oc 1 2 , which indicates much faster (1.5 times) expansion rate of the Universe. This clearly puts up severe 
problem in Large scale structure formation (LSS). 

1 2 

5. The solution in the early matter dominated era tracks t% power law, rather than standard t» , like any R~ n , 

with n > model, which again creates problem in structure formation and fitting WMAP data, as discussed 
earlier [8]. 

6. The parametric values which fit SNIa data, pushes the redshift at decoupling (zd ec ), which corresponds to 
3.75 x 10 -4 Billion years in the standard model, to very early epoch (zd ec ~ 200). However, if this data is fixed 
a-priori, then transition redshift (deceleration to acceleration) occurs at unacceptably large red-shift. 

Thus, unless supplemented by Einstein-Hilbert term, i?2 term alone suffers from serious disease. Since, Noether 
symmetry of F(R) theory of gravity does not admit anything other than R~s , in the Robertson- Walker minisu- 
perspace model. So, either we have to abandon Noether symmetry or try to involve more degrees of freedom in 
the theory, by incorporating either scalar field or working with some anisotropic minisuperspace model. This we 
pose in a future communication. 
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